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We show that equations governing pulsations of superfluid neutron stars (NSs) can be split into 
two sets of weakly coupled equations, one describing the superfluid modes and another one, the 
normal modes. The coupling parameter s is small, \s\ ~ 0.01 — 0.05, for realistic equations of state 
(EOSs). Already an approximation s = is sufficient to calculate the pulsation spectrum within 
the accuracy of a few percent. Our results indicate, in particular, that emission of gravitational 
waves from superfluid pulsation modes is suppressed in comparison to that from normal modes. 
The proposed approach allows to drastically simplify modeling of pulsations of superfluid NSs. 



PACS numbers: 97.60.Jd, 47.75. +f, 97.10.Sj, 47.37.+q 

Introduction. - - The pulsations of NSs can be ex- 
cited either due to internal instabilities or owing to ex- 
ternal perturbations. Currently the detectors that will 
be able, according to preliminary estimates, to register 
gravitational waves from pulsating NSs, are under de- 
velopment [l|. For the correct interpretation of future 
observations it is necessary to have a well developed the- 
ory of NS pulsations. The formulation of such a the- 
ory is complicated by the fact that at a temperature 
T < 10 8 - 10 10 K, baryons in the internal layers of NSs 
become superfluid. Thus, to model pulsations one has 
to employ superfluid hydrodynamics which is much more 
complicated than the ordinary one, describing "normal" 
(nonsuperfluid) matter. 

For the first time the global pulsations of superfluid 
NSs were analyzed by Lindblom and Mendell in 1994 0. 
Considering a simple model of a Newtonian star they nu- 
merically found two distinct classes of pulsation modes: 
(i) normal modes which practically coincide with the cor- 
responding modes of a normal star; and, (ii) superfluid 
modes in which the matter pulsates in such a way that 
the mass current density approximately vanishes. The 
subsequent numerical studies of various pulsation modes 
(the literature is vast; see, e.g., Refs. [3, 4] and references 
therein) confirmed the result of Ref. [2] though general 
explanation of this result has not yet been proposed [5| . 
In this work we give such an explanation. In addition, we 
present an approximate scheme which allows to greatly 
simplify calculations of pulsating superfluid NSs. In what 
follows the speed of light c = 1 . 

Superfluid hydrodynamics. - For simplicity, we 
consider NS cores composed of neutrons (n) , protons (p) , 
and electrons (e). We also assume that protons are nor- 
mal while neutrons are superfluid in some region of a 
NS core. As demonstrated in Ref. [fj possible admixture 
of other particle species (e.g., muons) and proton super- 
fluidity do not affect our principal results. Finally, we 
first consider a nonrotating NS. Effects of rotation will 
be incorporated later in the text. 

It is well known that in superfluid matter, several in- 
dependent motions with different velocities may coexist 
without dissipation [7]. In our case the system is fully 



described by two four- vectors, u M and w? n y The vec- 
tor m m is the velocity of electrons and protons as well as 
"normal" neutrons; the vector wf > arises from additional 
degrees of freedom associated with neutron superfluid- 
ity. In the nonrelativistic limit the spatial components 
of the four-vector w^ •, are related to superfluid veloc- 
ities V sn of the Landau-Khalatnikov theory Q by the 
equality Wr n ) = m n (V sn — u), where m n is the neutron 
mass; u is the spatial component of the "normal" four- 
velocity u^. The electron j(\ , proton j? ,, and neutron 

j? s current densities are expressed through the vectors 
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n e u^, jfa = n p U> J 
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Y nn w^ n y Here m is the number density of 



particles I = n, p, or e. The expression for j? ■. consists 
of two terms reflecting the fact that both normal and su- 
perfluid liquid components contribute to neutron current 
density. The coefficient Y nn has been calculated in Ref. 
[lfj |: it is a relativistic analogue of superfluid density of 
neutrons p sn . In the nonrelativistic limit Y nn — psn/m^. 
This coefficient depends on T and increases steadily from 
for normal matter (when T > T cn , where T cn is the neu- 
tron critical temperature) to n n /fJ, n f° r entirely super- 
fluid matter (T = 0). Here and below, \ii is the chemical 
potential for particles I = n, p, or e. 

To proceed further, we assume that: (i) the quasincu- 
trality condition holds both for equilibrium and pulsat- 
ing matter, n p = n e ; and, (ii) an unperturbed NS is 
in beta-equilibrium, i.e. the disbalance of chemical po- 
tentials Sjj, = ji n — jj,p — fi e = 0. It is convenient then 
to formulate the system of nondissipative hydrodynamic 
equations using the baryon number density rib = n n + n p 
and Sfi as independent variables Q . The system consists 
of (1) the continuity equations for baryons and electrons 
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where the baryon current density is j 
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njyU^ + Y nn W( y (2) Einstein equations 
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with the energy-momentum tensor T^ w — (P + e) u^u v + 
Pg^ + Y nn (wfaw{ n) + H n w<{ n) u v +fi n w , ( n) u^; (3) the 

potentiality condition for superfluid motion of neutrons 

dv [W(n)ii + VnUfj] - 9 p [W(n) v + /J, n Uv] = 0, (3) 

and (4) the second law of thermodynamics 
de = u„ dn h -Sudn P + T dS + Y nn d (w" n) W( n ) a ) /2. (4) 

In Eqs. ©-O, R^" ', R, and g^ are Ricci tensor, scalar 
curvature, and metric tensor, respectively; d^ = d/dx^; 
G is the gravitation constant; P = —e + \i n n b — 8[j,n e + 
T S is the pressure; e and S are the energy and en- 
tropy densities, respectively. All the thermodynamic 
quantities are defined in the comoving frame in which 
u'' = (1 ,0,0, 0). This imposes an additional constraint 



on w 
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Ufl w^ n) = 0. The solution to Eqs. ©-© in 
the superfluid region should be matched with that in the 
residual (normal) region of a star. Equations, describ- 
ing pulsations of normal matter can be obtained from 
the system HJ-Q if one puts Y nn = and ignores the 
condition §5§. 

Linear Approximation. — In this work we assume 
that in the unperturbed NS W( n \ — 0, i.e. velocities of 
normal and superfluid components coincide. For a nonro- 
tating NS this simply means that both components are at 



rest. Then it follows from the constraint 



\n) 



Othat 



w^\ vanishes in equilibrium, W/t-, — (because u° ^ 0). 
As a consequence, the terms in Eq. (|4|) and in the expres- 
sion for T 1 "' ', underlined with one line, are quadratically 
small. Similarly, the terms which depend on T and un- 
derlined twice, are small in the strongly degenerate mat- 
ter of NSs and can be omitted [9] . Because of the same 
reasons one can consider, e.g., the quantities P and 8/i as 
functions of only n b and n e (and neglect their dependence 
on scalars w^^w^n^ and T). 

Now we make use of the energy-momentum conser- 
vation law T.^f = which can be derived from Eq. ©. 
Composing a vanishing combination T.^, u +u^ u u T™ and 
subtracting from it Eq. ^ multiplied by n b « M , one gets, 
with the help of Eq. Q and the expression for P, 

—n e d^Sfi — u^ u v n e d^Sfi — n e Sfi u v {u^-v 

+{g^ v + u^Uu) u a X v a + X v Uu .„ + A M uY v 

-m u v [dvW^n ~ <9 M W(„)„] = 0. (5) 

Here A M = ii n Y nn w 1 ?, . The obtained "superfluid" equa- 
tion is very attractive because each term in it depends 
either on 8fi or itff ■>. Both these quantities are small in 
a slightly perturbed matter (and vanish in equilibrium). 
This means that in the linear approximation Eq. (J5)) does 
not depend explicitly on the perturbations of the metric 
g^v and the four- velocity u M . Thus, one can replace g^ u 
and u^ in Eq. ([5]) by their equilibrium values. For a non- 
rotating NS the spatial components of Eq. ([5]) can be 
rewritten in a remarkably simple form (J = 1, 2, 3) 



where we assumed that w?, depends on time t as w'f, ~ 

exp(iwi). In Eq. ([6]) all the quantities except for 8/j. and 
are taken in equilibrium. Near the equilibrium, the 



i r 



(n) _ 

function 8fi(n b , n e ) can be expanded in the Taylor series 
and presented, in the linear approximation, as 



Sfi = n e (d5(i/dn e ) (z Di + D 2 ). 



(7) 



Here z = [nbdSfi/dnb]/[n e ddfj,/dn e }; D\ = drib /rib; D 2 = 
5n e /n e . The symbol 5 in front of some quantity denotes a 
deviation of this quantity from its equilibrium value. The 
dimensionless functions D\ and D 2 can be found from Eq. 



(HJ) and depend on 



u 1 - 1 , and g^ u . Thus, generally, 



= ]? b) /n b = u» + Y 
instead of u^ and w?-, 
jjfj- = u m anc | \y" = q (since in equilibrium ivf 



Eq. ([5]) is not independent and should be solved together 
with Einstein equations ([2]) . In the linear approximation 
Eq. (J2J can be written in the following symbolic form: 
8{R tMV - 1/2 g» v R) = -8ttG 8T» V . The left-hand side of 
this equation contains only perturbations of metric. To 
write out the right-hand side it is convenient to introduce 
new independent variables, the four-velocity of baryons 

u» 

„ w , - 

Note that for an unperturbed star 

(») " °)' 
The same is also true for a pulsating NS if the matter 

is nonsuperfluid (because then Y nn = 0). Employing the 

new variables an expression for 8T^ V takes the form 

ST"" = {8P + 8e)U tl U" 

+(P + s)(U» 8U U + U v 8U^) + 6PgV + P8g^.(8) 

Here the quantities P, e, U^ and g^ v are taken at equi- 
librium. As follows from Eq. Q the variation 5e equals 
Se = fi n 5n b and depends on 8U^ and 8g^ v (and is inde- 
pendent of W). The variation 5P can be expanded in 
analogy with Eq. ([7]), 



SP = n b (dP/dn b )(D 1 + sD 2 ), 



(9) 



iU){Hn Ynn ~ n b ) W( n )j = n e dj(y/-g Sfl), (6) 



where the function D\ depends on 5U^ and 8g^ v , and D 2 
depends on the difference (SU^ — 8W^) and 8g^ v [see the 
definitions for D\ and D 2 , and Eq. (JTJ]. The parameter 
s, hereafter referred to as the "coupling parameter", is 
given by s = (n e dP/dn e )/(n b dP/dn b ). 

Superfluid and normal modes. - - If s = then 
Eq. © for ST^ V does not depend on 8W 11 and has exactly 
the same form as in the absence of superfluidity. In that 
case the Einstein equations (and boundary conditions to 
them) coincide, in the linear approximation, with the cor- 
responding equations for normal matter. They can be 
solved separately from the "superfluid" Eq. © so that 
the solution (the spectrum of eigenfrequencies uj and the 
eigenfunctions 8U^) will be indistinguishable from that 
for a nonsuperfluid star. 

Let us now focus on the following question. Assume 
that s still vanishes. Is it possible for a NS to oscil- 
late on a frequency which is not an eigenfrequency of a 
normal star? Suppose that it is indeed the case. Then 
the linearized Einstein equations will be satisfied only if 
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FIG. 1: The coupling parameter s versus density p for EOSs 
PAL HH, PA [H, APR pj|, SLy4 [l|, and LNS [3. EOSs 
of the PAL family differ by the symmetry energy (models I, 
II, or III) and by the value of the compression modulus, 180 
or 240 MeV. The models with the compression modulus 120 
MeV are not plotted since they give too small maximum NS 
masses that seem to contradict observations [16|| . Note that 
the recently measured mass M = (1.97 ± O.O4)M [l7J] of the 
millisecond pulsar PSR J1614-2230 further rules out PA EOS 
and all PAL EOSs except for PAL 1240 and PAL II240. 

6U 1 * = and Sg 1 *" = 0. As follows from Eq. (JTJ, in 
this case Srib — (i.e. D\ = 0), while D 2 depends only 
on 5W^ (or, in other words, on W^, since If = in 
equilibrium). In particular, for a nonrotating NS 



D 2 = [(d n e /n e ) W j + W ] A/{iujU°) 



(10) 



Here j — 1, 2, and 3; all the quantities, except for T4 7 - 7 , 
are taken in equilibrium; when calculating the covariant 
derivative one should use the unperturbed metric. Eqs. 
d?]) and (fTUj) allow to formulate Eq. © purely in terms 
of W^ . A boundary condition to this equation, W± = 0, 
also depends only on W J and can be obtained from the 
requirement that the baryon current density j^ is con- 
tinuous through the normal-superfhiid interface (W ± is 
the component of a vector W° , perpendicular to the in- 
terface). Thus, Eq. (j5]) is self-contained and can be solved 
independently of Eq. ^. Its solution (eigenfrequencies 
and eigenfunctions W° ) describes superfluid modes which 
were first considered in Ref. [2j and do not have an ana- 
logue for a normal star. To our best knowledge, the strik- 
ing properties of such modes have not been discussed for 
a realistic model of a general relativistic NS at finite T. 
First of all, the superfluid pulsation modes do not perturb 
metric (Sg^ — 0) and hence cannot emit gravitational 
waves. In addition, because for these modes <5L™ = and 
5rib = 0, the variations of jffi. and P vanish, 5j? b \ = 
and 6P = [see Eqs. (JXJ) and ©]. As a consequence, 
pulsations are localized entirely in the superfluid region 
of a star. In particular, they do not go to the NS surface. 



In the consideration above we supposed that s = 0. 
Yet, it is clear that superfluid and normal modes should 
remain approximately decoupled also at small but finite 
s. As follows from Fig. 1, s is indeed small for realistic 
EOSs and changes, on the average, from —0.01 to —0.05 
[18J . Taking into account that the parameter z in Eq. ([7]) 
is z ~ — 1 for the same EOSs, it is easy to show that for 
normal modes D\ > D 2 [then the second term in Eq. © 
is much smaller than the first one], while for superfluid 
modes D\ ~ s D 2 [then the first term in Eq. (JTJ) is much 
smaller than the second one]. Generally, the exact solu- 
tion of linearized pulsation equations JpJ and © can be 
presented as a series in parameter s [6|. However, since 
s is very small, the approximation of noninteracting Eqs. 
@ and ^ considered above (hereafter "zero approxi- 
mation" ) is already sufficient to calculate the pulsation 
spectrum within the accuracy ~ s (i.e., a few percent). 

Example: Radial pulsations. - Let us illustrate 
the obtained results with an example of a radially pulsat- 
ing NS with the mass M = IAMq. We consider a simple 
NS model which was analyzed in detail in Ref. [8J . In that 
paper it was assumed that the redshifted critical temper- 
ature of neutrons T^ is constant throughout the stellar 
core, T™ = 6 x 10 8 K. The results of approximate calcula- 
tion of pulsation spectrum are illustrated in Fig. 2a. The 
spectrum is calculated in zero approximation in s. In the 
figure, the pulsation frequency u (in units of ojq = c/i?NS, 
where i?NS — 12.17 km is the circumferential radius of a 
star) is plotted as a function of internal redshifted stellar 
temperature T°° for 3 normal (solid lines) and 6 super- 
fluid (dashes) pulsation modes. At T°° > T^ only the 
normal modes (I, II, and III) survive since then the star 
is nonsuperfluid. For comparison, in Fig. 2b we present 
the exact solution to the system of linearized equations 
© and ([S]). The first 6 modes are shown by alternate 
solid and dashed lines. The spectrum was not plotted in 
the shaded region. All other notations are the same as 
in Fig. 2a. It is easy to see that the structure of both 
spectra is very similar. However, there is one principal 
difference. Instead of crossings of superfluid and normal 
modes in Fig. 2a, we have avoided crossings in Fig. 2b. 
At these points the superfluid mode becomes normal and 
vice versa. Such avoided crossings are not described in 
approximate treatment (Fig. 2a) because when frequen- 
cies of superfluid and normal modes are close to each 
other, Eqs. ([2) and ^ become strongly interacting and 
cannot be considered as independent. For comparison, 
we plot both spectra in Fig. 2c. The exact solution is 
shown by solid lines, dashes correspond to the approxi- 
mate solution. Other notations are the same as in Figs. 
2a and 2b. On average, the approximate solution differs 
from the exact one by ~ 1.5 — 2%. For normal modes 
the difference becomes smaller with increasing of T . In 
this case the number of "superfluid" neutrons decreases 
(Y n „ -t 0), consequently, W J = Y nn w^ n) /n b -> and 
zero approximation works better and better. 

Taking into account rotation. — Rotation leads to 
formation of Feynman-Onsager vortices inside NSs with 
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FIG. 2: Frequency u in units of Wo versus T$° = T°°/(10 8 K) for various pulsation modes, (a) approximate spectrum; (b) 
exact spectrum; (c) approximate (dashed lines) and exact (solid lines) spectra. For more details see the text. 



the interspacing distance ~ 10~ 2 — 1CP 4 cm. The hydro- 
dynamic equations averaged over the volume containing 
large amount of vortices formally have the same form 
as in their absence [l9| (if we neglect the small con- 
tribution of vortices to the internal energy density of 
matter). The only exception is the potentiality condi- 
tion JjJJ) that should be replaced by u v {d^,[wr n ) v + fJ"n u v ] 
-9 y [w(„) M +/i„u M ]} = OfjtuW", where the tensor M „ 
is specified in Ref. [6| and is responsible for the inter- 
action between the normal and superfiuid component. 
It can be found from the requirement that the entropy 
does not decrease. Because of this condition the new 
term rib M „ W v appears in the right-hand side of Eq. 
([5]). Since this term depends on a small quantity W^, all 
our reasoning about decoupling of superfiuid and normal 
modes remain valid for rotating NSs as well. 

Comparison with previous works. - For com- 
parison we choose two papers, Refs. [20] and j4(, since at 
first sight it is not clear whether our results complement 
or contradict the conclusions drawn in these references. 

The authors of Ref. [2(| considered a model of Newto- 
nian star at T — 0. They demonstrated that superfiuid 
modes decouple from the normal modes only for an ide- 
alized case of nonstratified NSs, for which n e /rib = const 
throughout the stellar core. 

This result does not contradict ours because one can 
show that the neutron-star matter is nonstratified only if 
9P(7ib,n e )/9n e = (that is s — 0). As follows from our 
analysis, in the latter case superfiuid and normal modes 
are indeed strictly decoupled. 

The second conclusion made in Ref. [20| is based on 
the observation that for most of the neutron-star mod- 
els the stellar matter is stratified. Using this observa- 
tion the authors of Ref. [2fj] argued that generally there 
should be no clear distinction between the superfiuid and 
normal modes, or, in other words, equations describing 
superfiuid- and normal-type pulsations are strongly in- 
teracting. 

This conclusion is not correct because, as we demon- 
strated earlier in this work, the real coupling parameter 



s can be small even for strongly stratified NSs (and is 
indeed small for realistic EOSs). 

Now let us discuss the results of Ref. [J]. This paper 
analyzed gravitational radiation from superfiuid nonro- 
tating NSs at T — in the frame of the general relativity. 
It was argued that superfiuid modes must radiate gravi- 
tational waves in practically all situations, with intensity 
of radiation comparable to that from the normal modes 
(unless an EOS has a very specific form satisfying Eq. 
(74) of Ref. 0). 

When modeling the neutron-star pulsations the au- 
thors of Ref. [4j used toy-model EOSs that give com- 
pletely unrealistic values for the coupling parameter s. 
In particular, we found that their most realistic model II 
gives s ~ 0.1 at the center and s — oo at the superfluid- 
normal interface. Moreover, because their EOSs are ar- 
tificial, they were forced to relax an assumption of chem- 
ical equilibrium in the core. As it is demonstrated in the 
present paper, the latter assumption is very important for 
the decoupling of modes and cannot be ignored. Thus, 
it is not surprising that our results disagree with the re- 
sults of Ref. |4j ; when s is not small, superfiuid modes can 
be as effective in radiating gravitational waves as normal 
modes. 

In the end, it is worth mentioning one more result of 
Ref. 4]. In that paper it is claimed that any (nonradial) 
pulsation mode must emit gravitational waves unless an 
EOS satisfies some specific criterion [their Eq. (74)]. We 
checked that this criterion is not equivalent and does not 
follow from our criterion s = 0, which is a necessary 
condition for decoupling of superfiuid modes from metric. 

Conclusion. - Summarizing, equations describing 
pulsations of superfiuid NSs can be split into two sys- 
tems of weakly coupled equations. The coupling pa- 
rameter s of these systems is small for realistic EOSs, 
\s\ ~ 0.01 — 0.05. One system of equations describes 
normal modes, another one - superfiuid modes. Already 
zero approximation in parameter s (when the systems 
are fully decoupled) is sufficient to calculate the pulsa- 
tion spectrum with an accuracy of a few percent. In this 



approximation the normal modes coincide with ordinary 
modes of nonsuperfluid NS, while superfluid modes do 
not perturb metric, pressure, baryon current density and 
are localized in superfluid region of a star. Note that 
an emission of gravitational waves by superfluid modes 
is possible only in the next (first) order of perturbation 
theory in s. Thus, it should be suppressed in comparison 
to gravitational radiation from the normal modes. 

Our finding that superfluid modes do not appear at 
the NS surface and do not emit gravitational waves in 
the s — limit indicate that these modes should be very 
difficult to observe at small but finite s. This means that 
observational properties of a pulsating superfluid star and 
a normal star of the same mass should be very similar, 
so that it will be very hard to discriminate one from the 
other. 

The obtained results explain numerical calculations 



0, 3 and suggest simple perturbative (in parameter s) 
scheme which drastically simplifies the problem of cal- 
culation of the pulsation spectrum for superfluid NSs. 
The presented approach allows to easily take into account 
realistic EOSs, dissipation, various composition of mat- 
ter, temperature effects, baryon superfluidity, density- 
dependent profiles of critical temperatures, and rotation 
of NSs. In more detail these issues will be discussed else- 
where 6] . 

Acknowledgements. - We thank D.P. Barsukov, 
A.I. Chugunov, and D.G. Yakovlev for valuable com- 
ments. This research was supported by the Dynasty 
Foundation, Ministry of Education and Science of Rus- 
sian Federation (Contract No. 11.G34.31.000I with SPb- 
SPU and leading scientist G.G. Pavlov), RFBR (Grant 
No. 11-02-00253-a), and by FASI (Grant No. NSh- 
3769.2010.2). 



[1] N. Andersson, V. Ferrari, D. I. Jones et al., 

arXiv:0912.0384. [13] 

[2] L. Lindblom and G. Mendell, Astrophys. J. 421, 689 

(1994). [14] 

[3] U. Lee, Astron. Astrophys. 303, 515 (1995); L. Lind- 
blom and G. Mendell, Phys. Rev. D61, 104003 (2000); [15] 
S. Yoshida and U. Lee, Phys. Rev. D67, 124019 (2003); 
L.-M. Lin, N. Andersson, and G. L. Comer, Phys. Rev. [16] 
D78, 083008 (2008). 

[4] N. Andersson, G. L. Comer, and D. Langlois, Phys. Rev. [17] 
D66, 104002 (2002); 

[5] See, however, Refs. [201. |21j. where in two special cases [18] 
and under a number of simplified assumptions the au- 
thors exactly decouple, in the nonrelativistic limit, su- 
perfluid degrees of freedom from the normal ones. 

[6] M. E. Gusakov and E. M. Kantor, in preparation. 

[7] I. M. Khalatnikov, An Introduction to the Theory of Su- [19] 
perflmdity (Addison- Wesley, New York, 1989). 

[8] M. E. Gusakov and N. Andersson, Mon. Not. R. Astron. 
Soc. 372, 1776 (2006). 

[9] M. E. Gusakov, Phys. Rev. D76, 083001 (2007). [20] 

[10] M. E. Gusakov, E. M. Kantor, and P. Haensel, Phys. Rev. 

C79, 055806 (2009); Phys. Rev. C80, 015803 (2009). [21] 

[11] M. Prakash, J. M. Lattimer, and T. L. Ainsworth, Phys. 

Rev. Lett. 61, 2518 (1988). 
[12] D. Page and J. H. Applegate, Astrophys. J. 394, L17 



(1992). 

A. Akmal, V. R. Pandharipande, and D. G. Ravenhall, 
Phys. Rev. C58, 1804 (1998). 

E. Chabanat, P. Bonche, P. Haensel, J. Meyer, and R. 

Schaeffer, Nucl. Phys. A635, 231 (1998). 

L. G. Cao, U. Lombardo, C. W. Shen, and N. V. Giai, 

Phys. Rev. C73, 014313 (2006). 

J. M. Lattimer and M. Prakash, Phys. Rep. 442, 109 

(2007). 

P. B. Demorest, T. Pennucci, S. M. Ransom, M. S. E. 

Roberts, and J. W. T. Hessels, Nature, 467, 1081 (2010). 

The maximum value \s\ ~ 0.1 is reached for EOS PAL 

II180 at p ~ 10 15 g cm -3 (this region is not shown in Fig. 

1). At such high densities hyperons are likely to appear 

in NS matter so that the simple semiphenomenological 

EOS which allows only for nucleons, is not applicable. 

I. L. Bekarevich and I. M. Khalatnikov, Zh. Eksp. Teor. 

Fiz. 39, 1699 (1960) [Sov. Phys. JETP 12, 1187 (1961)]; 

G. Mendell and L. Lindblom, Ann. Phys. 205, 110 

(1991). 

R. Prix and M. Rieutord, Astron. Astrophys. 393, 949 

(2002); 

B. Haskell, N. Andersson, and A. Passamonti, Mon. Not. 
R. Astron. Soc. 397, 1464 (2009). 



